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ABSTRACT 

We identify high density regions of supercluster size in high resolution A-body simulations of a repre- 
sentative volume of three Cold Dark Matter Universes. By assuming that (1) the density and peculiar 
velocities of baryons trace those of the dark matter, and (2) the temperature of plasma is proportional 
to the velocity dispersion of the dark matter particles in regions where the crossing times is smaller 
than the supercluster free-fall time, we investigate how thermal motions of electrons in the intra-cluster 
medium and peculiar velocity of clusters can affect the secondary anisotropies in the cosmic microwave 
background (CMB). We show that the thermal effect dominates the kinematic effect and that the largest 
thermal decrements are associated with the most massive clusters in superclusters. Thus, searching for 
the presence of two or more close large CMB decrements represents a viable strategy for identifying 
superclusters at cosmological distances. Moreover, maps of the kinematic effect in superclusters are 
characterized by neighboring large peaks of opposite signs. These peaks can be as high as ~ 10^K at 
the arcminute angular resolution. Simultaneous pointed observations of superclusters in the millimeter 
and submillimeter bands with upcoming sensitive CMB experiments can separate between the thermal 
and kinematic effect contributions and constrain the evolution of the velocity field in large overdense 
regions. 

Subject headings: cosmic microwave background - large scale structure of the Universe - galaxies: 
clusters: general - gravitation - methods: A-body simulations 



1. INTRODUCTION 

Thermal motions of ionized gas in the intra-cluster 
medium and the bulk flow of clusters produce sec- 
ondary temperature fluctuations in the cosmi c microwave 
background (CM B) on the arcminute scale ( [gunyaev &j 
Zel'dovich 1970, 1972). For a given cluster velocity and 
mass, these thermal and kinematic effects (tSZ and kSZ ef- 
fects, hereafter) are independent of redshift. Typical line- 
of-sight peculiar velocities of galaxy clusters, estimated 
from distance indic ators to cluster galaxy members, are 
~ 300 km s" 1 (e.g. |Giovanelli et al. 1998| ). These veloc- 
ities can easily be ~ 50% larger for clu sters in superclus- 
ter environments flColberg et al. 1999 ). Therefore, CMB 
anisotropies caused by the kSZ effect in superclusters will 
be substantially enhanced. 

Ideally, one would like to estimate the magnitude of SZ 
anisotropies using hydro-dynamical simulations. Unfortu- 
nately, these simulations do not currently have the neces- 
sary dynamical range to study, at the same time, the large 
scale motions of clusters and the dynamics of the gas inside 
them. Therefore, we resort to high resolution A-body sim- 
ulations of dark matter particles only. By assuming that 
the motion and distribution of electron gas follow the dark 
matter, these simulations can be used to estimate the SZ 
anisotropies in the CMB. Here, we show that superclus- 
ters have a clear signature in both tSZ and kSZ maps: in 
tSZ maps, superclusters appear as regions of neighboring 
large CMB decrements; in kSZ maps, supercluster regions 



contain large neighboring peaks of opposite sign. The first 
property can provide an easy strategy for indentifying su- 
perclusters at any redshift; after removing the dominating 
tSZ effect contribution, the detectable kSZ peaks provide 
a direct measure of the velocity field in superclusters. 

2. METHOD 

To explore the non-linear evolution of the density and 
velocity field of superclusters, we use the A-body simula- 
tions of thr ee variants of the Cold Da rk Matter (CDM) 
model (see Kauffmann et al. 1999] for details): two 



low density models with a 141 J h J Mpc 3 comoving box 5 
(OCDM with n Q = 0.3; ACDM with n = 0.3, A = 0.7) 
and a high density model with a 85 3 /i -3 Mpc 3 comoving 
box (rCDM with f2o = 1). All models contain 256 3 parti- 
cles and are normalized to the present abundance of galaxy 
clusters. 

In order to locate superclusters in the simulation, we 
compute the density p(x) on a 256 3 grid from the particle 
positions at z = with the cloud-in-cell (CIC) interpo- 
lation scheme. We then smooth the density field with a 
gaussian window of width AhT x Mpc, compute the aver- 
age overdensity profiles S(< r) around local maxima and 
search for regions with S(r < 8/i _1 Mpc) = <5g > S s , where 
5 S is a parameter. The chosen window width identifies su- 
perclusters efficiently: smaller widths also identify isolated 
clusters; larger widths miss low density superclusters. Be- 
cause we set the same initial random phases in all our 



simulations, they yield the same large scale matter distri- 
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bution at z = 0. For example, in the OCDM model at 
z = 0, we find 45, 15, 7, 3, and 1 regions with S s = 2, 
3, 4, 5, and 6, respectively. The most dense region has 
<5g = 6.8. Hereafte r, we define as super clusters the regions 
with 5g > 4 (e.g. Einasto et al. 1997j ). These superclus- 
ters have crossing times larger than the Hubble time and 
they arc not in dynamical equilibrium. Clusters within 
the simulation box are identified with a friends-of-friends 
algorithm applied to the dark matter particles. 

The tSZ and kSZ effects depend on the temperature and 
line of sight velocity of the intra-cluster medium, and on 
the number density of free electrons. Therefore we have 
to relate these quantities to the positions and velocities 
of the dark matter particles in the simulations. We first 
compute the density and velocity fields in the simulation 
box on a 2048 3 grid with the CIC interpolation scheme. 
This grid preserves the 30ft-~ 1 kpc spatial resolution of our 
simulations in the high density regions we are interested 
in, namely superclusters, filaments, clusters and groups of 
galaxies. A coarser grid would artificially lower the magni- 
tude of the SZ fluctuations but not their relative location 
and sign. We then assume that at any point x within the 
cluster the plasma is thermalized: fcT e (x) = <7 2 (x)m p /2, 

where T e is the gas temperature, crjj the local one dimen- 
sional velocity dispersion of the dark matter particles, m p 
the proton mass, and k the Boltzmann constant. We also 
assume that the electron number density, n e (x), is related 
to the mass density p(x) by n c (x) = [p(x)/m p ]fib/fio, 
where fib is the background baryonic density parameters 
in terms of the critical density. All calcu lations in this Let- 
ter ar c done w ith fib = 0.0125/i -2 (e.g. Smith, Kawano & 
Malai fcy 1993| ). Note however, that clusters may contain 
more bary ons than the average baryon fracti on in the Uni- 
verse (e.g. Fukugita, Hogan fc Peebles 1998| ), so we might 
be underestimating the actual magnitude of the SZ effects 
by a factor of a few. 

Under these assumptions, the hot electron gas causes 
tSZ temperature decrements in the CMB spectrum 
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where x — hv /kT r , T r is the radiation temperature, L the 
simulation box proper size, ot the Thomson cross-section, 
m c the electron mass, and c the speed of light. Numerical 
values of the tSZ decrement presented here are computed 
in the Rayleigh- Jeans limit x — > 0. Line-of-sight peculiar 
velocities vu produce kSZ temperature fluctuations 
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With broad band observations, one can separate the ther- 
mal and kinematic effects and use the latter to measure 
the peculiar velocity of clusters (Sunyaev fc Zcl'dovicrj 



1980|: |Rcphaeli fc Lahav 199l| ; [Hachnclt fc Tegmark 1996 
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We obtain two dimensional maps by projecting the sim- 
ulation box along a random line of sight. Hereafter, we 
present results for all models projected along the same di- 
rection. We also compute maps of the Thomson optical 
depth tt = (fib/fio)°T Jq p(l) l 'm^dl , the line of sight ve- 



locity dispersion (er 2 ) = o?<(l) p(l)dl / p(l)dl, and the 
line of sight velocity field (v\\) — v\\ (l)p(l)dl/ p{l)dl. 

3. RESULTS 

Figure [j] shows the tSZ decrements in the OCDM box 
at z = 0. The largest decrements are in high density 
regions: superclusters appear as neighboring large decre- 
ments. This signature is characteristic of supercluster re- 
gions. Because the tSZ effect is independent of redshift, 
and, on arcminute angular scales, the tSZ decrements 
of the most massive clusters in a supercluster are larger 
than the cosmological fluctuations, neighboring large CMB 
decrements clearly indicate the presence of superclusters 
at cosmological distances. 

Figure shows the kSZ map corresponding to Figure 
There are many pairs of close fluctuations with oppo- 
site sign. The large scale structure clearly appears as a 
fuzzy filamentary network of small spots, which are the 
dark matter halos constituting the filaments. These are 
only a few Mpc long and are substantially curved. Be- 
cause of this curvature, filaments have relatively small op- 
tical depth and do not appreciably contribute to the CMB 
anisotropies with the kSZ effect. The largest kSZ contribu- 
tions comes from high density regions (Figure |||) . All the 
peaks in the tt map which show kSZ fluctuations larger 
than 5/iK are within regions with 5s > 2. There are 41 
peaks with tt > 4 x 10~ 3 of which 12 are outside regions 
with <5g > 2. There are no significant contributions to 
the kSZ fluctuations from peaks with tt < 1CU 3 . Peaks 
with t t > 1(T 5 obey the relation tSZ(^K) = -2 x 10 6 r|. 
Note that the bulk flow of the entire simulation box, when 
included, may increase the magnitude of the kSZ fluctua- 
tions. For example, in OCDM, the rms velocity of spheres 
of diameter lAOh^ 1 Mpc is ~ 90 km s _1 , and some kSZ 
fluctuations in Figure @ would increase by ~ 15%. 

Figure |] shows a supercluster region. Because the clus- 
ters are approaching each other, the kSZ peaks have oppo- 
site signs: this is the characteristic signature of the veloc- 
ity field within superclusters. The maps shown in Figure 
H] are 22.2 arcmin on a side and have resolution 0.4 ar- 
cmin; the two largest tSZ decrements have separation of 
5 arcmin. Examples of neighboring positive and negative 
kSZ peaks are also shown in Figure H for a supercluster at 
z = as it would appear in our three cosmological mod- 
els. The angular dimension of these maps depends on the 
supercluster distance. The tCDM, ACDM and OCDM 
models yield fluctuations with increasing magnitude. A 
full comparison of the magnitude of the SZ effects in our 
cosmological models will be discussed elsewhere. Here, we 
emphasize that, in all models, superclusters present neigh- 
boring tSZ decrements and kSZ fluctuations with opposite 
sign. Therefore our supercluster identification method is 
independent of the cosmology. 

Table 1 lists the properties of the two largest tSZ decre- 
ments in a few dense superclusters at different redshifts in 
the OCDM box. In some cases, superclusters have a pecu- 
liar feature: the kSZ fluctuations corresponding to the two 
tSZ decrements are not the largest fluctuations in the kSZ 
maps. We have searched for the kSZ fluctuations which are 
more than twice the kSZ fluctuation corresponding to the 
tSZ decrement. We find 9 such shifted kSZ fluctuations, 
out of the 3 superclusters at 5 different redshifts listed 
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Table 1 

Evolution of the Properties of Supercluster Regions 
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1 


0.0 


(9.4, 10.2) 


(1052,968) 


(138, -297) 


(-205. - 


193) 


(4,-10) 


0.5/0.1 


9.4/0.4 


0.6/0.3 


3.38 


1 


0.3 


(10.4,6.8) 


(1206, 1270) 


(447, -810) 


(-287, - 


217) 


(17. -20) 


0.7/0.4 


15.0/3.7 


1.1/0.7 


0.34 


1 


0.5 


(9.8, 10.6) 


(1201, 902) 


(361,449) 


(-284, - 


171) 


(12, 16) 


0.7/0.4 


12.3/3.4 


0.9/0.6 


0.72 


1 


0.7 


(14,6.5) 


(983, 1003) 


(-12, -159) 


(-277. - 


134) 


(0,-4) 


0.7/0.5 


10.0/4.1 


0.9/0.4 


4.28 


1 


1.1 


(13.2, 10.4) 


(744, 722) 


(106, 12) 


(-149, - 


110) 


(5,1) 


0.8/0.5 


6.1/2.3 


0.8/0.6 


3.45 


2 


0.0 


(8.7,3.9) 


(997,617) 


(-84, -254) 


(-172. - 


-29) 


(-3,-3) 


0.4/0.2 


5.9/1.4 


0.4/0.3 


1.48 


2 


0.3 


(10.1, 5.4) 


(948,551) 


(-11, -55) 


(-185, - 


-32) 


(o,-i) 


0.5/0.3 


6.2/0.8 


0.5/0.5 


4.48 


2 


0.5 


(10.9,8) 


(969,722) 


(-577, 174) 


(-206, - 


-85) 


(-23,5) 


0.5/0.4 


6.8/1.4 


0.8/0.5 


0.90 


2 


0.7 


(9.8,9.3) 


(745,753) 


(-515.9) 


(-111,- 


106) 


(-17,0) 


0.5/0.3 


4.6/0.7 


0.7/0.4 


2.20 


2 


1.1 


(11.6,9.5) 


(749,678) 


(-338, -15) 


(-131, - 


-89) 


(-12,0) 


0.6/0.4 


4.9/0.7 


0.6/0.4 


3.27 


5 


0.0 


(6.3,5.2) 


(789,718) 


(-367, 568) 


(-80. - 


55) 


(-8,10) 


0.3/0.1 


2.8/0.1 


0.4/0.1 


1.36 


5 


0.3 


(6.9,7.3) 


(742, 706) 


(-350,408) 


(-77. - 


75) 


(-8,10) 


0.4/0.1 


3.1/0.1 


0.4/0.1 


2.03 


5 


0.5 


(8.4,7.8) 


(710, 703) 


(430, -346) 


(-86. - 


78) 


(12,-9) 


0.4/0.2 


3.3/1.0 


0.5/0.2 


1.82 


5 


0.7 


(8.8,6.9) 


(719,613) 


(462, -327) 


(-91,- 


53) 


(14,-8) 


0.5/0.3 


3.4/1.5 


0.6/0.2 


1.45 


5 


1.1 


(10.9, 7.3) 


(666,579) 


(427, -109) 


(-101, - 


-49) 


(17,-3) 


0.6/0.3 


3.4/1.0 


0.7/0.4 


4.20 



Evolution of the properties of 8 X 8h Mpc supercluster regions centered on the supercluster positions at z — 0. Values are from maps of the OCDM 
model with a comoving spatial resolution of 0.2/i _1 Mpc. Column (1) supercluster identification number in order of decreasing ovcrdensity; column (2) 

rcdshift; columns (3)-(7) optical depth, velocity dispersion, peculiar velocity, tSZ and kSZ fluctuations of the largest (left entries in brackets) tSZ 
decrement and second largest (right entries) tSZ decrement; columns (8)— (10) rms values including/excluding two squares of 2h~~ l Mpc on a side each 
centered on one of the two largest tSZ decrements; (11) projected comoving separation between the two largest tSZ decrements. 



in Table 1. For example, in supercluster #2 at z = 0, 
there are two kSZ fluctuations of — 7/iK and 6/iK at a pro- 
jected distance ~ l/i -1 Mpc and ~ 3/i _1 Mpc from the 
tSZ decrement — 172^tK, whose corresponding kSZ fluctu- 
ation is — 3/iK. These shifted kSZ peaks originate from 
clusters which are different from the clusters yielding the 
two largest tSZ decrements and are not necessarily relaxed. 

The presence of large kSZ fluctuations and tSZ decre- 
ments that do not overlap is confirmed by the rms of the 
map fluctuations within the supercluster regions. In Ta- 
ble 1 we list two rms's, including or excluding two squares 
of 2h^ 1 Mpc on a side, each centered on one of the two 
largest tSZ decrements. The tSZ rms drops substantially 
when these two squares are excluded. However, this is 
not always the case for the tt and kSZ rms's. Therefore, 
within superclusters, there may be halos with large pecu- 
liar velocities and sufficiently high Thomson optical depth 
which cause kSZ fluctuations larger than the fluctuations 
due to the most massive clusters. We finally note that in 
our two low density models, amplitudes of both the tSZ 
and kSZ effects do not show any substantial evolution in 
the redshift range < z <; 0.7. At higher redshifts, how- 
ever, the average amplitudes decrease appreciably. 

4. CONCLUSION 

The recent development in the microwave detector tech- 
nology, the planned MAP, Planck Surveyor, and FIRST 



experiments, large ground arrays (e.g. ALMA and in 
Antarctica), and various balloon projects will provide high 
precision measurements of CMB temperature fluctuations 
on small scales. Here, we have seen that peculiar motions 
of clusters in supercluster environment produce detectable 
arcminute anisotropies. These kSZ fluctuations come in 
the form of peaks of opposite signs but they are an order 
of magniture smaller than temperature decrements esti- 
mated from the tSZ effect in the Rayleigh- Jeans limit. 

In single-band low frequency CMB observations, super- 
clusters at any redshift can be identified as well sepa- 
rated areas of neighboring non-gaussian large decrements 
superimposed on cosmological and foregrounds contribu- 
tions which are more uniformly distributed. Follow-up ob- 
servations in the millimeter and submillimcter bands can 
disentangle the thermal and kinematic contributions. We 
find that the largest thermal and kinematic fluctuations 
can originate from different clusters within the same su- 
percluster and, therefore, might not always be perfectly 
overlapping. 

The GIF iV-body simulations were carried out at the 
Computer Center of the Max-Planck Society in Garching 
and at the EPPC in Edinburgh, as part of the Virgo Con- 
sortium project. A. N. is grateful to the MPA for the 
hospitality. 
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Fig. 1. — Map of the tSZ decrements (black spots) in the OCDM box at z = 0. The comoving resolution is 0.2h~ Mpc. To increase the 
contrast, decrements larger than 3 times the rms= 1.2/iK are saturated. The 15 most dense regions are located by squares of 8h~ 1 Mpc on 
a side. All but one of the high density region pairs, whose squares overlap on the map, are nearby in 3D: the separations between the two 
maxima defining each region are < 12/i -1 Mpc. The pair 2-8 only is a projection effect, with a separation between the maxima of Alh~ 
Mpc. 



Fig. 2. — Map of the kSZ fluctuations corresponding to Figure pi White (black) spots correspond to positive (negative) fluctuations. To 
increase the contrast, fluctuations larger than 3 times the rms= 0.12/^K are saturated. The 7 most dense superclusters are located by squares 
of S/i" 1 Mpc on a side. Superclusters # 1 and # 3, and superclusters # 6 and # 5 are nearby in 3D: the separations between the two maxima 
defining each supercluster are 9 and 8 Mpc, respectively. 




Fig. 3. — Amplitude of the tSZ and kSZ fluctuations corresponding to peaks in the OCDM tt map at z = 0. Upper panels show peaks 
outside the 45 regions with <5g > 2. Lower panels show peaks within these 45 regions. The dashed lines in the left panels show that the tight 
tSZ-TT relation is well described by the law tSZ(/^K) = — 2 X 10 6 t|,. 



Fig. 4. — OCDM tt, tSZ and kSZ maps in supercluster # 5 at z = 0.5. Maps have comoving size 8xSh 2 Mpc 2 and 0.2h 1 Mpc comoving 
resolution. 



Fig. 5. — kSZ maps in the supercluster region # 4 in our three cosmological models at z = 0. These models have the same initial random 
phases, so they yield the same large scale matter distribution at z = 0. Maps have comoving size 8 X 8h~ 2 Mpc 2 and 0.2h _1 Mpc comoving 
resolution. 
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